Understanding the impact of mobility on COVID-19 spread: A hybrid gravity-metapopulation model of COVID-19.

The outbreak of the severe acute respiratory syndrome coronavirus 2 started in Wuhan, China, towards the end of 2019 and spread worldwide. The rapid spread of the disease can be attributed to many factors including its high infectiousness and the high rate of human mobility around the world. Although travel/movement restrictions and other non-pharmaceutical interventions aimed at controlling the disease spread were put in place during the early stages of the pandemic, these interventions did not stop COVID-19 spread. To better understand the impact of human mobility on the spread of COVID-19 between regions, we propose a hybrid gravity-metapopulation model of COVID-19. Our modeling framework has the flexibility of determining mobility between regions based on the distances between the regions or using data from mobile devices. In addition, our model explicitly incorporates time-dependent human mobility into the disease transmission rate, and has the potential to incorporate other factors that affect disease transmission such as facemasks, physical distancing, contact rates, etc. An important feature of this modeling framework is its ability to independently assess the contribution of each factor to disease transmission. Using a Bayesian hierarchical modeling framework, we calibrate our model to the weekly reported cases of COVID-19 in thirteen local health areas in Metro Vancouver, British Columbia (BC), Canada, from July 2020 to January 2021. We consider two main scenarios in our model calibration: using a fixed distance matrix and time-dependent weekly mobility matrices. We found that the distance matrix provides a better fit to the data, whilst the mobility matrices have the ability to explain the variance in transmission between regions. This result shows that the mobility data provides more information in terms of disease transmission than the distances between the regions.

A generalized distributed delay model of COVID-19: An endemic model with immunity waning.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been spreading worldwide for over two years, with millions of reported cases and deaths. The deployment of mathematical modeling in the fight against COVID-19 has recorded tremendous success. However, most of these models target the epidemic phase of the disease. The development of safe and effective vaccines against SARS-CoV-2 brought hope of safe reopening of schools and businesses and return to pre-COVID normalcy, until mutant strains like the Delta and Omicron variants, which are more infectious, emerged. A few months into the pandemic, reports of the possibility of both vaccine- and infection-induced immunity waning emerged, thereby indicating that COVID-19 may be with us for longer than earlier thought. As a result, to better understand the dynamics of COVID-19, it is essential to study the disease with an endemic model. In this regard, we developed and analyzed an endemic model of COVID-19 that incorporates the waning of both vaccine- and infection-induced immunities using distributed delay equations. Our modeling framework assumes that the waning of both immunities occurs gradually over time at the population level. We derived a nonlinear ODE system from the distributed delay model and showed that the model could exhibit either a forward or backward bifurcation depending on the immunity waning rates. Having a backward bifurcation implies that $ R_c < 1 $ is not sufficient to guarantee disease eradication, and that the immunity waning rates are critical factors in eradicating COVID-19. Our numerical simulations show that vaccinating a high percentage of the population with a safe and moderately effective vaccine could help in eradicating COVID-19.

Association between close interpersonal contact and vaccine hesitancy: Findings from a population-based survey in Canada

Background Vaccine hesitancy threatens efforts to bring the coronavirus disease 2019 (COVID-19) pandemic to an end. Given that social or interpersonal contact is an important driver for COVID-19 transmission, understanding the relationship between contact rates and vaccine hesitancy may help identify appropriate targets for strategic intervention. The purpose of this study was to assess the association between interpersonal contact and COVID-19 vaccine hesitancy among a sample of unvaccinated adults in the Canadian province of British Columbia (BC). Methods Unvaccinated individuals participating in the BC COVID-19 Population Mixing Patterns Survey (BC-Mix) were asked to indicate their level of agreement to the statement, “I plan to get the COVID-19 vaccine.” Multivariable multinomial logistic regression was used to assess the association between self-reported interpersonal contact and vaccine hesitancy, adjusting for age, sex, ethnicity, educational attainment, occupation, household size and region of residence. All analyses incorporated survey sampling weights based on age, sex, geography, and ethnicity. Results Results were based on survey responses collected between March 8, 2021 and December 6, 2021, by a total of 4,515 adults aged 18 years and older. Overall, 56.7% of respondents reported that they were willing to get the COVID-19 vaccine, 27.0% were unwilling and 16.3% were undecided. We found a dose-response association between interpersonal contact and vaccine hesitancy. Compared to individuals in the lowest quartile (least contact), those in the fourth quartile (highest contact), third quartile and second quartile groups were more likely to be vaccine hesitant, with adjusted odd ratios (aORs) of 2.85 (95% CI: 2.02, 4.00), 1.91(95% CI: 1.38, 2.64), 1.78 (95% CI: 1.13, 2.82), respectively. Conclusion Study findings show that among unvaccinated people in BC, vaccine hesitancy is greater among those who have high contact rates, and hence potentially at higher risk of acquiring and transmitting infection. This may also impact future uptake of booster doses.

Assessing the potential impact of immunity waning on the dynamics of COVID-19 in South Africa: an endemic model of COVID-19.

We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of backward bifurcation and bi-stability, where a stable disease-free equilibrium coexists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19.

The basic reproduction number of COVID-19 across Africa.

The pandemic of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) took the world by surprise. Following the first outbreak of COVID-19 in December 2019, several models have been developed to study and understand its transmission dynamics. Although the spread of COVID-19 is being slowed down by vaccination and other interventions, there is still a need to have a clear understanding of the evolution of the pandemic across countries, states and communities. To this end, there is a need to have a clearer picture of the initial spread of the disease in different regions. In this project, we used a simple SEIR model and a Bayesian inference framework to estimate the basic reproduction number of COVID-19 across Africa. Our estimates vary between 1.98 (Sudan) and 9.66 (Mauritius), with a median of 3.67 (90% CrI: 3.31-4.12). The estimates provided in this paper will help to inform COVID-19 modeling in the respective countries/regions.

Social Contacts and Transmission of COVID-19 in British Columbia, Canada.

Background: Close-contact rates are thought to be a driving force behind the transmission of many infectious respiratory diseases. Yet, contact rates and their relation to transmission and the impact of control measures, are seldom quantified. We quantify the response of contact rates, reported cases and transmission of COVID-19, to public health contact-restriction orders, and examine the associations among these three variables in the province of British Columbia, Canada. Methods: We derived time series data for contact rates, daily cases and transmission of COVID-19 from a social contacts survey, reported case counts and by fitting a transmission model to reported cases, respectively. We used segmented regression to investigate impacts of public health orders; Pearson correlation to determine associations between contact rates and transmission; and vector autoregressive modeling to quantify lagged associations between contacts rates, daily cases, and transmission. Results: Declines in contact rates and transmission occurred concurrently with the announcement of public health orders, whereas declines in cases showed a reporting delay of about 2 weeks. Contact rates were a significant driver of COVID-19 and explained roughly 19 and 20% of the variation in new cases and transmission, respectively. Interestingly, increases in COVID-19 transmission and cases were followed by reduced contact rates: overall, daily cases explained about 10% of the variation in subsequent contact rates. Conclusion: We showed that close-contact rates were a significant time-series driver of transmission and ultimately of reported cases of COVID-19 in British Columbia, Canada and that they varied in response to public health orders. Our results also suggest possible behavioral feedback, by which increased reported cases lead to reduced subsequent contact rates. Our findings help to explain and validate the commonly assumed, but rarely measured, response of close contact rates to public health guidelines and their impact on the dynamics of infectious diseases.

Effect of Human Mobility on the Spatial Spread of Airborne Diseases: An Epidemic Model with Indirect Transmission.

We extended a class of coupled PDE-ODE models for studying the spatial spread of airborne diseases by incorporating human mobility. Human populations are modeled with patches, and a Lagrangian perspective is used to keep track of individuals' places of residence. The movement of pathogens in the air is modeled with linear diffusion and coupled to the SIR dynamics of each human population through an integral of the density of pathogens around the population patches. In the limit of fast diffusion pathogens, the method of matched asymptotic analysis is used to reduce the coupled PDE-ODE model to a nonlinear system of ODEs for the average density of pathogens in the air. The reduced system of ODEs is used to derive the basic reproduction number and the final size relation for the model. Numerical simulations of the full PDE-ODE model and the reduced system of ODEs are used to assess the impact of human mobility, together with the diffusion of pathogens on the dynamics of the disease. Results from the two models are consistent and show that human mobility significantly affects disease dynamics. In addition, we show that an increase in the diffusion rate of pathogen leads to a lower epidemic.

Mathematical modeling of COVID-19 in British Columbia: An age-structured model with time-dependent contact rates.

Following the emergence of COVID-19 at the end of 2019, several mathematical models have been developed to study the transmission dynamics of this disease. Many of these models assume homogeneous mixing in the underlying population. However, contact rates and mixing patterns can vary dramatically among individuals depending on their age and activity level. Variation in contact rates among age groups and over time can significantly impact how well a model captures observed trends. To properly model the age-dependent dynamics of COVID-19 and understand the impacts of interventions, it is essential to consider heterogeneity arising from contact rates and mixing patterns. We developed an age-structured model that incorporates time-varying contact rates and population mixing computed from the ongoing BC Mix COVID-19 survey to study transmission dynamics of COVID-19 in British Columbia (BC), Canada. Using a Bayesian inference framework, we fit four versions of our model to weekly reported cases of COVID-19 in BC, with each version allowing different assumptions of contact rates. We show that in addition to incorporating age-specific contact rates and mixing patterns, time-dependent (weekly) contact rates are needed to adequately capture the observed transmission dynamics of COVID-19. Our approach provides a framework for explicitly including empirical contact rates in a transmission model, which removes the need to otherwise model the impact of many non-pharmaceutical interventions. Further, this approach allows projection of future cases based on clear assumptions of age-specific contact rates, as opposed to less tractable assumptions regarding transmission rates.

Assessing the potential impact of immunity waning on the dynamics of COVID-19 in South Africa: an endemic model of COVID-19

We developed an endemic model of COVID-19 to assess the impact of vaccination and immunity waning on the dynamics of the disease. Our model exhibits the phenomenon of backward bifurcation and bi-stability, where a stable disease-free equilibrium coexists with a stable endemic equilibrium. The epidemiological implication of this is that the control reproduction number being less than unity is no longer sufficient to guarantee disease eradication. We showed that this phenomenon could be eliminated by either increasing the vaccine efficacy or by reducing the disease transmission rate (adhering to non-pharmaceutical interventions). Furthermore, we numerically investigated the impacts of vaccination and waning of both vaccine-induced immunity and post-recovery immunity on the disease dynamics. Our simulation results show that the waning of vaccine-induced immunity has more effect on the disease dynamics relative to post-recovery immunity waning and suggests that more emphasis should be on reducing the waning of vaccine-induced immunity to eradicate COVID-19.

Quantifying transmissibility of SARS-CoV-2 and impact of intervention within long-term healthcare facilities.

Estimates of the basic reproduction number (R 0) for COVID-19 are particularly variable in the context of transmission within locations such as long-term healthcare (LTHC) facilities. We sought to characterize the heterogeneity of R 0 across known outbreaks within these facilities. We used a unique comprehensive dataset of all outbreaks that occurred within LTHC facilities in British Columbia, Canada as of 21 September 2020. We estimated R 0 in 18 LTHC outbreaks with a novel Bayesian hierarchical dynamic model of susceptible, exposed, infected and recovered individuals, incorporating heterogeneity of R 0 between facilities. We further compared these estimates to those obtained with standard methods that use the exponential growth rate and maximum likelihood. The total size of outbreaks varied dramatically, with range of attack rates 2%-86%. The Bayesian analysis provided an overall estimate of R 0 = 2.51 (90% credible interval 0.47-9.0), with individual facility estimates ranging between 0.56 and 9.17. Uncertainty in these estimates was more constrained than standard methods, particularly for smaller outbreaks informed by the population-level model. We further estimated that intervention led to 61% (52%-69%) of all potential cases being averted within the LTHC facilities, or 75% (68%-79%) when using a model with multi-level intervention effect. Understanding of transmission risks and impact of intervention are essential in planning during the ongoing global pandemic, particularly in high-risk environments such as LTHC facilities.

Assessing the impact of adherence to Non-pharmaceutical interventions and indirect transmission on the dynamics of COVID-19: a mathematical modelling study.

Adherence to public health policies such as the non-pharmaceutical interventions implemented against COVID-19 plays a major role in reducing infections and controlling the spread of the diseases. In addition, understanding the transmission dynamics of the disease is also important in order to make and implement efficient public health policies. In this paper, we developed an SEIR-type compartmental model to assess the impact of adherence to COVID-19 non-pharmaceutical interventions and indirect transmission on the dynamics of the disease. Our model considers both direct and indirect transmission routes and stratifies the population into two groups: those that adhere to COVID-19 non-pharmaceutical interventions (NPIs) and those that do not adhere to the NPIs. We compute the control reproduction number and the final epidemic size relation for our model and study the effect of different parameters of the model on these quantities. Our results show that there is a significant benefit in adhering to the COVID-19 NPIs.

How much leeway is there to relax COVID-19 control measures?

Following successful non-pharmaceutical interventions (NPI) aiming to control COVID-19, many jurisdictions reopened their economies and borders. As little immunity had developed in most populations, re-establishing higher contact carried substantial risks, and therefore many locations began to see resurgence in COVID-19 cases. We present a Bayesian method to estimate the leeway to reopen, or alternatively the strength of change required to re-establish COVID-19 control, in a range of jurisdictions experiencing different COVID-19 epidemics. We estimated the timing and strength of initial control measures such as widespread distancing and compared the leeway jurisdictions had to reopen immediately after NPI measures to later estimates of leeway. Finally, we quantified risks associated with reopening and the likely burden of new cases due to introductions from other jurisdictions. We found widely varying leeway to reopen. After initial NPI measures took effect, some jurisdictions had substantial leeway (e.g., Japan, New Zealand, Germany) with > 0.99 probability that contact rates were below 80% of the threshold for epidemic growth. Others had little leeway (e.g., the United Kingdom, Washington State) and some had none (e.g., Sweden, California). For most such regions, increases in contact rate of 1.5-2 fold would have had high (> 0.7) probability of exceeding past peak sizes. Most jurisdictions experienced June-August trajectories consistent with our projections of contact rate increases of 1-2-fold. Under such relaxation scenarios for some regions, we projected up to ∼100 additional cases if just one case were imported per week over six weeks, even between jurisdictions with comparable COVID-19 risk. We provide an R package covidseir to enable jurisdictions to estimate leeway and forecast cases under different future contact patterns. Estimates of leeway can establish a quantitative basis for decisions about reopening. We recommend a cautious approach to reopening economies and borders, coupled with strong monitoring for changes in transmission.

Importance of COVID-19 vaccine efficacy in older age groups.

IMPORTANCE: An effective vaccine against SARS-CoV-2 will reduce morbidity and mortality and allow substantial relaxation of physical distancing policies. However, the ability of a vaccine to prevent infection or disease depends critically on protecting older individuals, who are at highest risk of severe disease. OBJECTIVE: We quantitatively estimated the relative benefits of COVID-19 vaccines, in terms of preventing infection and death, with a particular focus on effectiveness in elderly people. DESIGN: We applied compartmental mathematical modelling to determine the relative effects of vaccines that block infection and onward transmission, and those that prevent severe disease. We assumed that vaccines showing high efficacy in adults would be deployed, and examined the effects of lower vaccine efficacy among the elderly population. SETTING AND PARTICIPANTS: Our mathematical model was calibrated to simulate the course of an epidemic among the entire population of British Columbia, Canada. Within our model, the population was structured by age and levels of contact. MAIN OUTCOME(S) AND MEASURE(S): We assessed the effectiveness of possible vaccines in terms of the predicted number of infections within the entire population, and deaths among people aged 65 years and over. RESULTS: In order to reduce the overall rate of infections in the population, high rates of deployment to all age groups will be critical. However, to substantially reduce mortality among people aged 65 years and over, a vaccine must directly protect a high proportion of people in that group. CONCLUSIONS AND RELEVANCE: Effective vaccines deployed to a large fraction of the population are projected to substantially reduce infection in an otherwise susceptible population. However, even if transmission were blocked highly effectively by vaccination of children and younger adults, overall mortality would not be substantially reduced unless the vaccine is also directly protective in elderly people. We strongly recommend: (i) the inclusion of people aged 65 years and over in future trials of COVID-19 vaccine candidates; (ii) careful monitoring of vaccine efficacy in older age groups following vaccination.

Quantifying the impact of COVID-19 control measures using a Bayesian model of physical distancing.

Extensive non-pharmaceutical and physical distancing measures are currently the primary interventions against coronavirus disease 2019 (COVID-19) worldwide. It is therefore urgent to estimate the impact such measures are having. We introduce a Bayesian epidemiological model in which a proportion of individuals are willing and able to participate in distancing, with the timing of distancing measures informed by survey data on attitudes to distancing and COVID-19. We fit our model to reported COVID-19 cases in British Columbia (BC), Canada, and five other jurisdictions, using an observation model that accounts for both underestimation and the delay between symptom onset and reporting. We estimated the impact that physical distancing (social distancing) has had on the contact rate and examined the projected impact of relaxing distancing measures. We found that, as of April 11 2020, distancing had a strong impact in BC, consistent with declines in reported cases and in hospitalization and intensive care unit numbers; individuals practising physical distancing experienced approximately 0.22 (0.11-0.34 90% CI [credible interval]) of their normal contact rate. The threshold above which prevalence was expected to grow was 0.55. We define the "contact ratio" to be the ratio of the estimated contact rate to the threshold rate at which cases are expected to grow; we estimated this contact ratio to be 0.40 (0.19-0.60) in BC. We developed an R package 'covidseir' to make our model available, and used it to quantify the impact of distancing in five additional jurisdictions. As of May 7, 2020, we estimated that New Zealand was well below its threshold value (contact ratio of 0.22 [0.11-0.34]), New York (0.60 [0.43-0.74]), Washington (0.84 [0.79-0.90]) and Florida (0.86 [0.76-0.96]) were progressively closer to theirs yet still below, but California (1.15 [1.07-1.23]) was above its threshold overall, with cases still rising. Accordingly, we found that BC, New Zealand, and New York may have had more room to relax distancing measures than the other jurisdictions, though this would need to be done cautiously and with total case volumes in mind. Our projections indicate that intermittent distancing measures-if sufficiently strong and robustly followed-could control COVID-19 transmission. This approach provides a useful tool for jurisdictions to monitor and assess current levels of distancing relative to their threshold, which will continue to be essential through subsequent waves of this pandemic.